Exact Casimir-Polder potentials: 
interaction of an atom with a conductor-patched dielectric surface 
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We study the interaction between a neutral atom or molecule and a conductor-patched dielectric 
surface. We model this system by a perfectly reflecting disc lying atop of a non-dispersive dielectric 
half-space, both interacting with the neutral atom or molecule. We assume the interaction to be 
non-retarded and at zero temperature. We find an exact solution to this problem. In addition 
we generate a number of other useful results. For the case of no substrate we obtain the exact 
formula for the van der Waals interaction energy of an atom near a perfectly conducting disc. We 
show that the Casimir-Polder force acting on an atom that is polarized in the direction normal to 
the surface of the disc displays intricate behaviour. This part of our results is directly relevant 
to recent matter-wave experiments in which cold molecules are scattered by a radially symmetric 
object in order to study diffraction patterns and the so-called Poisson spot. Furthermore, we give 
an exact expression for the non-retarded limit of the Casimir-Polder interaction between an atom 
and a perfectly-conducting bowl. 

PACS numbers: 37.30.-K, 41.20.Cv, 42.50.Pq 
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I. INTRODUCTION 

The importance of the Casimir-Polder interaction in 
cold- atom physics is by now a well established fact and 
the theory of atom-surface interactions is well developed. 
There are theoretical tools available, both numerical [l[ 
and analytical @, Q , for calculating dispersion forces be- 
tween atoms and polarizable objects for arbitrary geome- 
tries and a variety of optical properties of the surfaces. 
However, practical calculations of Casimir-Polder forces 
tend to be complicated because, in one way or another, 
they require solutions to complex electromagnetic scat- 
tering problems. In order to make calculations manage- 
able and render final results to be of practical use, sim- 
plifying assumptions of some sorts need to be introduced. 
One possibility is to focus attention on the so-called non- 
retarded regime of the Casimir-Polder interaction, also 
called the van der Waals regime. In this regime one 
assumes the speed of light to be infinite and therefore 
any scattering problems whose solutions are required for 
the computation of the Casimir-Polder force become es- 
sentially electrostatic problems. This regime is applica- 
ble to atoms or molecules whose distance to the surface 
is much smaller than the wavelength of their dominant 
dipole transition. Although going to the non-retarded 
regime somewhat restricts the applicability of the final 
results, the benefits are enormous; electrostatics is con- 
siderably simpler than electrodynamics plus there are 
numerous scenarios where the non-retarded calculations 
are adequate and provide accurate estimates of Casimir- 
Polder forces. For example, the spacings between molec- 
ular energy levels are relatively close, so that the wave- 
lengths of dominant dipole transitions in molecules tend 
to be relatively large. Thus, the interaction between cold 
molecules and surfaces often falls into the non-retarded 
regime. Since the mathematical framework of electro- 
statics is very well investigated, one is able to obtain 



exact analytic solutions for some non-trivial geometries 
that bear more resemblance to realistic experimental set- 
ups than the oversimplified geometries often studied in 
the literature. 

In this paper we respond to the growing need of ex- 
perimental physicists for simple and easy-to-use expres- 
sions enabling efficient calculations of the Casimir-Polder 
forces in realistic geometries that tend to arise in modern 
experiments. The question that we think requires urgent 
attention is this: what is the interaction of an atom with 
a dielectric surface onto which some conducting struc- 
tures have been deposited? This is a very common ex- 
perimental scenario e.g. with atom chips. Inspired by 
recent progress in graphene technologies Q one might 
even imagine an atom interacting with graphene flake or 
ribbon supported by a dielectric substrate. Such struc- 
tures are now routinely made and it is just a matter of 
time until they find their way into cold-atom physics. 

With this motivation in mind, we model a conductor- 
patched dielectric surface interacting with a neutral atom 
or molecule by a perfectly reflecting disc supported by 
a non-dispersive dielectric half-space interacting with a 
neutral localized quantum system, represented by a fluc- 
tuating electric dipole moment. We assume the interac- 
tion to be non-retarded and at zero temperature. 

In the following we shall find an exact solution to 
the problem of a perfectly reflecting disc atop a non- 
dispersive dielectric, which we think is a quite remark- 
able result and to the best of our knowledge the first of 
its kind in the literature. In addition we generate a num- 
ber of other useful results. Most notably we obtain an 
exact formula for the van der Waals interaction energy of 
an atom near a perfectly conducting disc. We anticipate 
that this result will prove useful for current matter-wave 
experiments where cold molecules are scattered on radi- 
ally symmetric objects in order to study the diffraction 
patterns and in particular the so-called Poisson spot, a 
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phenomenon well-known from standard optics. We also 
derive an expression for the Casimir-Poldcr force acting 
on an atom close to the edge of a conducting half-plane 
deposited on a dielectric substrate. This formula can be 
applied to the case of an atom near the edge of a con- 
ducting ribbon. Furthermore, in the appendix we give 
an exact expression for the Casimir-Polder interaction 
between an atom and a perfectly-conducting bowl. 



II. NON-RETARDED CASIMIR-POLDER 
POTENTIAL 

At close range the Casimir-Poldcr interaction between 
a neutral atom and a perfectly reflecting material struc- 
ture is non-retarded and can be worked out by methods of 
electrostatics. Taking the atom to be a point-like electric 
dipole located at ro, one can show that the interaction 
energy with a nearby surface can in rectangular coordi- 
nates be written as p 
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where G//(r, r') is the homogeneous part of the Green's 
function, to be explained further below. The sum runs 
over the three components of the electric dipole moment 
operator fa. The complete Green's function G(r, r') is 
the electrostatic potential at r due to a unit point charge 
at r' and satisfies 



V 2 G(r,r') = <S (3) (r-r') 



(2) 



with appropriate boundary conditions. For a perfect re- 
flector, the boundary conditions on G(r, r') are for it to 
vanish for any r on the surface of the reflector. For dielec- 
tric surfaces Maxwell's equations imply continuity con- 
ditions on the gradient of G(r,r'). The component of 
VG(r, r') that is parallel to the surface of the dielectric 
is required to be continuous across the interface, whereas 
the component normal to the surface of the dielectric is 
continuous when multiplied by the respective dielectric 
constant e(r) on either side of the interface. The homo- 
geneous part of the potential is the difference between 
the Green's function G(r,r') and its free-space equiva- 
lent, i.e. the potential of the same point charge in free 
space, 



G ff (r,r') =G(r,r')- 
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So far the problem is entirely classical and does not in- 
volve quantum electrodynamics. The quantum proper- 
ties of the atom are accounted for in the expectation 
values of the electric dipole moment operator (pf) = 
(j\pi\j), where |j) denotes the state of the atom, not 
necessarily its ground state. The difficulty of determin- 
ing the energy shift ([T]) lies in calculating the Green's 
function G(r,r') for the geometry of interest. 



III. CONDUCTING HALF-PLANE ON A 
DIELECTRIC SUBSTRATE 

A. Green's function 

In this section we derive the Green's function of the 
Poisson equation for a semi-infinite conducting half-plane 
lying on top of a dielectric substrate, as depicted in Fig. 
[T] We work in cylindrical coordinates r = (p, <fi, z) and 
place the origin of the coordinate system on the edge of 
the conducting half-plane which is described by <j) = 0. 
The plane </> = ir describes the surface of the polarizable 
body that is not coated by the conductor. We require the 
electrostatic potential G(r, r') to vanish on the surface of 
the conducting half-plane, that is 



G(r,r') =0 for0 = {O,27r}. 



(4) 



Across the half-plane described by <fr — tt wc impose stan- 
dard dielectric continuity conditions, 
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G(r,r'), e((/))— G(r, r') continuous. 



(5) 

To find the solution of Eq. (|2|) that satisfies the bound- 
ary and continuity conditions described above, we work 
with the eigenfunctions ^ n (p, <j>, z) of the Laplace opera- 
tor in cylindrical coordinates and construct the Green's 
function in terms of an eigenfunction expansion (cf. e.g. 
Sect. Ill of Ref. (H). In the region < < 7r eigenfunc- 
tions that are regular for small p and that vanish on the 
half-plane <f> — are 
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FIG. 1. (Color online) A point charge at r' = (p ,<j>' ,z') near 
a perfectly reflecting half-plane lying on top of a dielectric 
half-space with constant and frequency independent dielectric 
function e = n 2 generates the potential given in Eq. (|f 0|) . Also 
shown are lines of constant potential. 
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In the region 7r < < 2n those that are regular for small 
p and vanish on the half-plane = 2ir are 



e lKZ J™ (kp)—= sm 2?r - 

2ir 2 V ^ V 2 



(7) 

The parameters /c, fc, and m in Eqs. ^ and ([7]) would 
not a priori need to be the same, but the requirement of 
continuity of the p and z derivatives, according to Eq. ([5]). 
forces them to be the same. In order to satisfy the simul- 
taneous continuity condition on e((b)d^ n (p,(f>, z)/d(f> one 
needs to have either 
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B. Energy shift 
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In the first case, m must therefore be an even integer, 
and in the second an odd integer. Putting the parts be- 
low and above — n together with the appropriate rel- 
ative factors so as to satisfy the continuity conditions at 
= 7r and then normalizing each cigenfunction, one can 
proceed to assembling the Green's function from these 
eigenfunctions and obtains in the region < < 7r and 
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and in the region 7r < <p < 2tt and < <// < 7r 
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Carrying out the integrations, as explained in detail in 
Sect. IV of Ref. Q , one then finds for the Green's function 
in the case of the source being above the material, that 
is for 0' e (0,7r), 
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n = 1.5 



FIG. 2. (Color online) Direction of the Casimir-Polder force 
acting on an neutral atom with isotropic polarizability placed 
near a conducting half-plane lying on a dielectric substrate 
with refractive index n = 1.5 (blue arrows). The approxi- 
mately horizontal lines are contours of constant van der Waals 
energy. 



With the Green's function (fT0|) determined, the calcu- 
lation of the energy shift is now straightforward. We sub- 
tract the free-space Green's function from Eq. fIU|) and 
plug the result into formula ([I]). For an atom located at 
(p, </), z), cf. Fig. [U we find that the non-retarded energy 
shift in the atom may be expressed as 

AE = - 2n 2 (J+l)e 1 3 >W> + E ^> + S '<^ 

(13) 

with the abbreviations 

5 cos0 [(n 2 + 1)tt - 20] [1 + sin 2 0] 



3 16p 3 sin 2 32p 3 sin 3 

1 cos \{n 2 + 1)tt - 20] [1 + cos 2 , 



48p 3 8p 3 sin 2 32p 3 sin 3 , 

1 cos0 [(n 2 + 1)tt - 20] 



24p 3 ' 16p 3 sin 2 . 



32p 3 sin 3 « 



In this form the expression for the energy shift facili- 
tates a clear-cut comparison with the result obtained 
in Ref. @ for the case of the atom interacting with 
a half-plane alone i.e. without any substrate present. 
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However, for practical purposes it is much more con- 
venient to express the shift in terms Cartesian coordi- 
nates where the directions of the unit vectors, in terms 
of which the dipole matrix elements (/if) are expressed, 
are position- independent. We set up the Cartesian co- 
ordinate system in such a way that the z direction is 
perpendicular to the surface and the x axis runs along 
the edge of the conducting sheet. Mathematically we set 
pcoscj) = y, psin0 — z, z = x in Eq. (|10p and calculate 
the energy shift again by using (Q}. We obtain 

[E x {^) + ^l) + ^M)] 



AE 



2ir 2 (n 2 + l)e 



with the abbreviations 
1 
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(17) 



For an atom with isotropic polarizability, 
which is the most common case, we get 

Z X (HI) + Zy(l4) + Z z (f4) = S iso (^ 2 ) With 

i y 



8(y 2 + z 2 ) 3/2 4z 2 (y 2 + z 2 ) 
n 2 n — 2 arctan (y/z) 



(18) 



which is a remarkably simple end result. It is straight- 
forward visualize the result by plotting the direction of 
the Casimir-Poldcr force, which is done in Fig. [2l As 
one would expect from physical intuition, there is a lat- 
eral component of the Casimir-Polder force due to the 
presence of the conducting coating. One can easily con- 
vince oneself, analytically or numerically, that the lateral 
component of the Casimir-Poldcr force pulling the atom 
towards the edge of the half-plane is dominated by the 
normal component of the force even for relatively low 
values of the index of refraction n. 



IV. KELVIN INVERSION 

The Kelvin inversion is a non-linear coordinate 
transformation, a reflection of space in a sphere of ra- 
dius S and centred at s; it is defined as follows: 



T\r] 



s 2 



r (r - s) + s. 



(19) 



Ref. [7j gives a very clear overview of the geometrical 
properties of this transformation. It is of interest here be- 
cause it preserves solutions of boundary-value problems 
of the Poisson equation, as we shall explain further below. 
The Green's function of the Poisson equation G(r, r') is a 
function of two variables: the observation point r and the 
source point r'. To apply the transformation (fl"§)) to the 
Green's function G(r, r') one applies it to both of its ar- 
guments. The so transformed Green's function is, when 
multiplied by an appropriate pre-factor, a solution to the 
Poisson equation in the transformed geometry, because 
one has [all] 



S 2 



r — s r — s 



r G(T[r],TM) 



<5 (3) (r-r'). (20) 



Therefore the transformation 
S 2 



( 15 ) G(r,r') 



r — s r — s 



G(T[r],T[r']) = G(r,r') (21) 



generates a new Green's function of the Poisson equa- 
tion G(r,r') in a new geometry from an already known 
Green's function G(r, r') in a geometry that is related 
to the new one by a Kelvin transformation. It generates 
a new solution to a boundary-value problem in a new 
geometry from the known solution in the original geom- 
etry, because, if the Green's function G(r, r') vanishes 
on some surface a, then the transformed Green's func- 
tion G(r,r') vanishes on the transformed surface a. In 
this way Green's functions for perfect reflectors of vari- 
ous shapes can be obtained from known Green's functions 
by various Kelvin transformations, i.e. by adjusting the 
position and the radius of the inversion sphere [l(J □J. 

For the present purposes it is important to note that, 
in general, the Kelvin transformation does not preserve 
dielectric boundary conditions, which is why it is nor- 
mally thought to be useful only for problems involving 
just perfect reflectors. The most obvious example of a 
case in which the Kelvin inversion does not preserve the 
solution of the Poisson equation is that of a dielectric 
sphere, in contrast to the potential of a point charge near 
a perfectly conducting sphere which may be obtained by 
applying the Kelvin inversion to the potential near the 
perfectly reflecting flat mirror. The potential near a flat 
mirror, perfectly reflecting or dielectric, is a sum of a 
free-space potential and its image in the mirror. The po- 
tential near a perfectly reflecting sphere also has the same 
structure, and the image charge and its location can be 
obtained from the flat-mirror solution by a Kelvin trans- 
formation. However, the potential of a point charge near 
a dielectric sphere cannot be written down as a sum of 
a free-space potential plus an image potential from just 
a single point image charge. The exact solution in this 
geometry is more complicated and involves a point im- 
age charge and an additional continuous line of image 
charges inside the sphere [l2j ■ Therefore the solutions in 
those two geometries, for the dielectric half-space and for 
the dielectric sphere, cannot be connected by the Kelvin 
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transformation (jT9|) . It is easy to understand why this 
is the case: a dielectric surface enforces different conti- 
nuity conditions on the gradients of the potential normal 
and parallel to it, but in general Kelvin transformations 
distort geometries and hence do not preserve normal and 
parallel directions. However, there are some special cir- 
cumstances under which the Kelvin inversion does in fact 
preserve normal and parallel directions at a surface, for 
example, for a cone whose tip is the centre of the inver- 
sion. This transformation is trivial as it maps a dielectric 
cone into itself but, as we shall see in the next section, it 
is nevertheless very useful as it can be applied to cones 
with conducting coatings. Note that a plane is a cone 
with 180° opening angle, and thus the geometry of Fig. 
for which the Green's function is given by Eq. (fT0|) . is 
a special case of such a cone. 



CONDUCTING DISC ON A DIELECTRIC 
SURFACE 



Green's function 
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FIG. 3. (Color online) A conducting disc lying on top of 
a dielectric substrate, which is a dielectric half-space with 
dielectric constant e — n 2 . 

In this section we are going to derive the Green's func- 
tion of the Poisson equation for the geometry of the con- 
ducting disc lying on top of a dielectric substrate, cf. Fig 
[31 We take the perfectly reflecting disc of diameter d to 
lie in the z = plane, centred at the origin. The elec- 
trostatic potential is required to vanish on the surface of 
the disc, i.e. 




FIG. 4. (Color online) Green's function (|10|l applies to the 
geometry of a conducting half-plane on a dielectric substrate 
where the perfect reflector occupies the y — PI x > 
plane (left). Application of a Kelvin inversion (|19|) with 
s = (— d, 0, 0) and S = d yields the Green's function for a 
conducting disc on a dielectric substrate (right). 



the Green's function with the above boundary and conti- 
nuity conditions by standard methods, such as an eigen- 
function expansion, would be a formidable, if not im- 
possible task. Instead, we are going to derive the re- 
quired Green's function by means of a Kelvin inversion, 
as described in the previous section, starting from the 
Green's function ([TO]) obtained in Sec. IIII Al for a half- 
space covered with a conducting half-plane. As explained 
at the end of Sec. IIV1 a Kelvin transformation preserves 
the correct continuity conditions on the surface of the 
dielectric provided the centre of the inversion sphere is 
placed somewhere on the surface of the dielectric sub- 
strate. We place the centre of the inversion sphere at 
s = (— d, 0, 0) and set its radius S = d. For this particular 
choice of s and S the transformation leaves the dielectric 
substrate unchanged while mapping the conducting half- 
plane a = {r e M 3 : y = n x > 0}, cf. Fig. H(lcft), into 
a disk ff = {rel 3 :i/ = 0n (x + d/2) 2 + z 2 < (d/2) 2 }, 
cf. Fig. [4] (right). This ensures the correct boundary and 
continuity conditions across the whole surface y = 0. 

To simplify the notation, we wish to center the coordi- 
nate system on the disc and orient the z-axis perpendic- 
ular to the sheet, as shown in Fig. [3l Thus we shift and 
rotate the axes according to x — > y — d/2, y — > z, z — > x. 
We are interested in the disc-geometry Green's function 
for the case z,z' > 0, so that we need the Green's func- 
tion (JTOj) for the case 4>, (j) 1 € (0, it). In that case Eq. ([12]) 
reads 



G(r, r') = for z = n x 2 + y 2 < 



4 



(22) 



where x, y are the coordinates in the plane of the disc. In 
the plane z = outside the disc we require that G(r, r') 
satisfies dielectric boundary conditions, that is, we re- 
quire that 

^ G( r, , | G (ry), e(z)§-G(r,ri) 

are continuous for z = n x 2 + y 2 > d 2 /A. Obtaining 



V- = !) V+ = sg n [sin (0 + <j>')} 



(23) 



It is advantageous to further rewrite rj + in terms of 
sin(0 + (j)') = sin((/)) cos(0') + cos(0) sin(</>') as this facil- 
itates the straightforward application of the Kelvin in- 
version. The procedure described above then yields the 
transformed Green's function in the same form as in Eq. 
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lUDI) but with 



with the abbreviations 



± d 2 zz' + 



P ~2 



z 2 + [p+- 



z' 2 + \p' - 



2 -i 



z' 2 + [p' + 




D T = vV + p' 2 - 2pp'cos(0 - </>') + (2 T z') 2 , (25) 
77+ = sgn[z(d 2 /4 - p' 2 - z' 2 ) + z'(d 2 /4 - p 2 - z 2 )], (26) 

and, of course, 77 = 1 . We emphasize that this result, 
i.e. Eq. (JlO]) with the insertions (|M] ) -([26 ]) . is valid for 
z, z' > 0, when the source and observation points are 
both in the empty space above the substrate, cf. Fig [3] 
The Green's function for the case of the source in vacuum 
but the observation point inside the material can also be 
obtained quite easily by considering rj± in Eq. (|T2|) and 
making appropriate amendments. 



B. Energy shift 
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FIG. 5. (Color online) Direction of the Casimir-Polder force 
acting on a neutral atom with isotropic polarizability placed 
near a conducting disc lying atop a dielectric substrate of 
refractive index n = 1.5. The approximately horizontal lines 
are contours of constant van der Waals energy. 



Knowing the Green's function we can now calculate 
the energy shift. We subtract the free-space potential 
from the result for the Green's function that we obtained 
in the Sec. IV Bl and substitute it in into formula ([1]). For 
an atom located at (p, <j>, z), cf. Fig. [3J we find that the 
non-retarded energy shift may be expressed as 
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where we have defined R± = yj{p ± d/2) 2 + z 2 . The 
result ([2T)I together with Eqs. ([50)1 is, as far as elec- 

trostatics is concerned, exact. It applies to atoms whose 
distance from the surface is much smaller than the wave- 
length of their dominant dipole transition, so that effects 
of retardation are unimportant. The simplicity of the an- 
alytic expressions in Eqs. (J28J) (J30J) makes it easy to eval- 
uate and plot the non-retarded Casimir-Polder force felt 
by an atom due to the presence of a dielectric substrate 
with a conducting circular patch. In Fig. [5] we show the 
direction of the force acting on an atom with isotropic 
polarizability. As expected there is a lateral component 
in the Casimir-Polder force that drags the atom towards 
the conducting patch. 

Furthermore, it is interesting to consider an atom po- 
larized in the z direction interacting with just a disc, 
i.e. without any substrate. To this end we take {p 2 } = 
(pA) = in Eq. (|27| . so that the whole of the energy 
shift comes from the function S z , and take the limit of 
no substrate, n — > 1. For such an atom that is polarized 
in the z direction the energy shift vanishes in the plane 
of the disc, and similarly to the interaction of an atom 
with a conducting sheet with a circular hole, as studied 
in Ref. [l]|, we expect some intricate behaviour of the 
Casimir-Polder force in this case. Since our results are 
exact we can easily visualize the force in such a situation, 
which is what we do in Fig. [Sj We note that along the 
edge of the disc there is a region where the component 
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FIG. 6. (Color online) Direction of the Casimir-Polder force 
acting on a neutral atom placed near a conducting disc. The 
atom is polarized in the direction perpendicular to the surface 
of the disc. 

of the Casimir-Polder force normal to the surface of the 
disc points away from the disc. This unusual behaviour of 
the Casimir-Polder force may potentially have an impact 
on matter-wave experiments which are sensitive to van 
der Waals interactions. For example, Ref. [l4| studies 
the interference patterns (Poisson spot) from molecules 
scattered by a disc. The theoretical calculations used in 
conjunction with such matter-wave experiments tend to 
rely on rather oversimplified ex pre ssions for the atom- 
obstacle dispersive interactions [15(. As we can see in 
Fig. [SI the Casimir-Polder interaction for a polarized 
molecule near a disc exhibits a complicated pattern, so 
that theoretical predictions using simplistic models may 
potentially differ appreciably from the true values that 
can be calculated on the basis of the exact results de- 
rived here, especially for beams of polarized molecules. 
Although Eqs. p8|) - ([30)l may seem complicated at first 
glance they are expressed in terms of only elementary 
functions and can thus be evaluated numerically without 
much effort. 



VI. CONCLUSIONS 

We have fulfilled our aim of calculating the exact van 
der Waals potential for a neutral atom interacting with a 
conducting disc supported by a dielectric half-space. To 
obtain the exact solution for this complicated but prac- 
tically very relevant geometry we first solved the Pois- 
son equation for the geometry of a conducting half-plane 
that is supported by a dielectric half-space, cf. Sec lIII Al 
Then, in Sec. IV Al we applied a Kelvin inversion to ob- 
tain the Green's function for a conducting circular disc 
on a dielectric substrate. The result for this Green's 
function is exact and, to the best of our knowledge, the 
first of its kind in the literature. It can be used to cal- 
culate static electric fields near dielectric surfaces with 
conducting coatings or patches, cf. Fig [7J The main 
goal of this paper was to determine the energy-level shift 
(van der Waals energy) for an atom in close proximity 
of such a structure. The remarkably simple result for 
the shift given by Eq. fl27|) and Eqs. pg j) -(|gn ]) allows 



one to estimate the Casimir-Polder interaction for atoms 
near dielectric substrates with conducting structures on 
top. One could imagine an electroplated and partially 
etched surface or a graphene flake on a Si02 substrate or 
a graphene quantum dot. In the limit of no substrate our 
result reduces to that for a flat disc. Detailed analysis of 
the atom-disc interaction reveals the intricate behaviour 
of the Casimir-Polder force acting on atoms or molecules 
polarized in the direction normal to the surface of the 
disc, cf. Fig. [SI This, as we pointed out at the end of 
Sec. IV B( may be of significance for matter-wave inter- 
ferometry experiments. 




FIG. 7. (Color online) Contours of constant electrostatic po- 
tential <j)(r) for a point charge above a dielectric half-space 
with a perfectly conducting disc on top. The index of refrac- 
tion of the substrate is n = 2. 

In Appendix [Al wc obtain the exact expression for the 
interaction of a neutral atom with a conducting spheri- 
cal bowl. Our interest in this geometry was ignited by 
nanocups, i.e. dielectric spheres that are half-coated in 
gold, which have been recently shown to enable second 
harmonic generation (l6j . However, our result does not 
include the dielectric filling of the nanocup and is there- 
fore applicable only approximately for weak dielectric 
substrates. 



Appendix A: Conducting spherical shell. Green's 
function and the energy shift 

The Green's function for a conducting disc without 
a substrate facilitates the derivation of yet another in- 
teresting exact solution in electrostatics, namely the po- 
tential due to a point charge in the presence of a con- 
ducting spherical bowl. In order to derive the Green's 
function in this geometry we apply a Kelvin inversion to 
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the Green's function for a conducting disc, which is ob- 
tained by taking the limit n — > 1 in the results of Sec. 
IV Al We place the disc of radius d in the z = plane 
and choose s = (0,0, d) as the centre of the Kelvin in- 
version and S = d as the radius of the inversion sphere, 
cf. Eq. (|19p and Fig. [5] We dispense with the details 





diz y 



FIG. 8. (Color online) Applying a suitable Kelvin inversion 
(|19[1 to the Green's function for a disc generates the Green's 
function for a conducting spherical bowl. Note that dielectric 
continuity conditions are not preserved by this transforma- 
tion, so that the procedure works only in the absence of the 
substrate, i.e. in the limit n — > 1. 



of the calculation and explicit formulae for the Green's 
function, which have already been reported in Ref. fiol ]. 
Here we are interested only in giving the energy-level 
shift for an atom interacting with a conducting spherical 
bowl. We take the centre of the bowl to be at the centre 
of the Cartesian coordinate system so that its surface is 
described by 



4 
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V z- 
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(A4) 

\\ ts ~\/ Z J 
■ y -+- z" — d /4, and the result for 5 y is ob- 
tained by the replacement x O y, that is, Sj,(x, y, z, d) = 
S 2; (j/, x, z, d). It may at first sight seem awkward to work 
in Cartesian coordinates in this geometry, but this in fact 
avoids mathematical difficulties which arise at the origin 
of the spherical coordinate system where the azimuthal 
angle is not well defined. Besides, it is useful to have 
dipole matrix elements expressed in terms of unit vectors 
whose orientation does not depend on their position. We 
plot the direction of the Casimir-Polder force acting on 
an atom with isotropic polarizability in Fig. [5J 

\ \ WW W1MIM11 WW ////// 




Then the energy shift can be expressed in Cartesian co- 
ordinates as 




FIG. 9. (Color online) Direction of the Casimir-Polder force 
acting on an atom with isotropic polarizability placed near a 
conducting spherical bowl (blue arrows). The black circular 
lines are contours of constant van-der-Waals energy, and the 
red semicircle is the conducting bowl. 



AE = rn^To [Ex{ ^ } + Sy< ^ } + Ez{fX ^ (A2) 
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